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ABSTRACT 

We investigate the effects of ambipolar diffusion and the Hall effect on the stability of weakly- 
ionized, magnetized planar shear flows. Employing a local approach similar to the shearing- 
sheet approximation, we solve for the evolution of linear perturbations in both streamwise- 
symmetric and non-streamwise-symmetric geometries using WKB techniques and/or numer- 
ical methods. We find that instability arises from the combination of shear and non-ideal 
magnetohydrodynamic processes, and is a result of the ability of these processes to influence 
the free energy path between the perturbations and the shear. They turn what would be sim- 
ple linear-in-time growth due to current and vortex stretching from shear into exponentially- 
growing instabilities. Our results aid in understanding previous work on the behaviour of 
weakly-ionized accretion discs. In particular, the recent finding that the Hall effect and am- 
bipolar diffusion destabilize both positive and negative angular velocity gradients acquires a 
natural explanation in the more general context of this paper. We construct a simple toy model 
for these instabilities based upon transformation operators (shears, rotations, and projections) 
that captures both their qualitative and, in certain cases, exact quantitative behaviour. 

Key words: instabilities - MHD - ISM: magnetic fields - ISM: jets and outflows - accretion, 
accretion discs 



1 INTRODUCTION 

The im portance of understand ing t he phys i cs of shear flows has been appreciated for well over a century, starting with the pioneering 
work of von Helmholtz 1 1868h and KelvinI 1871 ). The stability of these flows in a wide variety of situations has been thoroughly stud- 



ied and reviewed, most notably in Chandrasekhar's (1961) classic text. Since that time, observa tions have revealed that shear flows are 
commonplace in astro physical systems, playing likely role s in the stability an d coUimation of je ts ( e^.. Ferrari. Trussoni & Zaninettilll981 



objects (lPringle.1 1989,: | Bachiller„ 1996.) . rotating stars dCowling 



Balbus & Hawlev 



1951 



" 1^ Is J " 

Fiedler & Jones 1984; Begelman, Blandford & Reesll 19841). the solar corona je.g.. K arpen et al. 1993), bipolar out flows from young stellar 



Goldreich & Schuberli,1967.) . accretion discs i Papaloizou & Linlll995 



19981) . and even the Earth's magnetopause ' McKenzie 1 19701) . In nearly all of these systems, magnetic fields play a vital 



role in determining the structure and evolution. There are, however, systems where the importance of magnetic fields remain unclear, due to 
poor ionization and th erefore weak coupling betwee n the dominant ion and ne utral species and the m agnetic field. These include molecular 
cloud s and their cores ^Mouschovias & Ciolekll 19991). galacti c molecular discs i Blaes & Balbus|[l994 ). protostellar accretion discs i Gammie 
19961 : iHawlev & S tone" 1998: Fromang . Terguem & Balbus"2002'). protostellar outflows and disc winds (Wardle & Konigl 1993), dwarf 
nova disks i Gammie & Menou 1998 : Sano & Storiq , 2003). shock waves in dense molecular clouds tWardle 1991: Draine & McKe&T993l: 
Roberge & Ciolekl2007l) . and protoplanetarv discs jSano & Mivamal 1999l : ISano et al.l200d : ISalmeron & Wardlel2005l : lchiang & Murrav-Clav 
2007|) . In these systems, various non-ideal magnetohydrodynamic (MHD) effects may come into play. For example, the neutral particles may 



drift relative to the ions in a process referred to as ambipolar diffusion. If, on the other hand, the ions drift relative to the electrons (and 
thereby the magnetic field), then the Hall effect occurs. 

Given the high occurrence of shear flows and low levels of ionization in various astrophysical environments, it is not surprising that 
examples can be found when the two coincide. Perhaps the most obvious exa mples are outflows from star-forming regions, where shear 
instabilities can occur at the interface between the jet and the ambient material dWatson et al ] l2004h . The resulting turbulent boundary layer 



* Email: mkunz@uiuc.edu 



2 M.W. Kunz 



can transfer linear momentum from the jet to the ambient medium, suggesting a possible mechanism for entrainment. Even at the launching 
point of these outflows, there may be differing velocity profiles in the ion and neutral fluids (see fig. 2 of Wardle & Konigl 1993), conditions 
ripe for both ambipolar diffusion and shear. Another notable example concerns the interface between the magnetically-active and -inactive 
('dead') regions in protoplanetary discs. It has been suggested that the well-coupled outer layers of protoplanetary discs will be subject t o 
the magnetorotational instability (MRI; Balbus & Hawley 1991), while the shielded midplane layers will remain dormant i Gammie 19961) . 
In this situation, angular momentum can be effe ctively transported in the outer layers but not in the midplane, setting up a velocity profile in 
the vertical (away from the midplane) direction i Fleming & Stonell2003 : Fromang & Nelsonlbood) . 

Here we undertake a study of a magnetized, planar shear layer in the presence of ambipolar diffusion and the Hall effect, and fin d 
that instability arise s when shear and non-ideal MHD effects act in concert. The instabilities are similar to those found bv lWardle Jl999l) , 
Balbus & Terguem 1 2001 ), Kunz & BalbusI 1 2004 ), and DeschI 1 2004 ) in the context of the MRI. Here we show that they are actually much 



more general. Our approach not only provides us with a clearer physical picture of the instabilities, free from the complications of rotation, but 
also brings qualitatively new results. In the case of ambipolar diffusion, instability arises from the combination of shear and the anisotropic 
nature of its wave damping. When the Hall effect is present, epicyclic-like motions set up by electromagnetic (whistler) waves couple to the 
background shear. When the handedness of these waves is opposite to that of the shear, instability occurs. It is notable that n either of these 



proce sses relies on rotational k inematics. In fact, these instabi lities are neither versions of the Kelvin-Helmholtz instability dWatson et al 



2004 ), nor versions of the MRI i Balbus & Hawlevll991 , 1 9921)1) . despite their reliance on the presence of both shear and magnetic fields. The 



unstable modes are a result of the ability of non-ideal MHD processes to open new pathways for the fluid to tap into the free energy of shear. 

An outline of the paper is as follows. In i)2.1l we discuss the formulation of the problem and give the basic equations to be solved. We 
then c onsider the evolution of Eulerian perturbations to these equations in comoving, local Lagrangian coordinates jGoldreich & Lvnden-Bell 
19651) . After deriving two coupled equations for the evolution of the relevant magnetic field eigenvectors in ^2.21 we solve them both 
analytically and numerically in several different situations. In ijs] we first restrict our attention to the effects of ambipolar diffusion on 
the stability of the system. After a detailed analysis and discussion, we then examine the effects of Hall electromotive forces (HEMFs) in 
Section|4] Contact with prior work is emphasized. In ijs] we construct and analyse a toy model that captures all the salient features of these 
instabilities. Section|6]summarizes our findings and conclusions. 



2 FORMULATION OF THE PROBLEM 



2.1 Basic Equations 

The equations describing a non-ideal MHD system, in the limit of negligible ion and electron inertia, are the continuity equation, 

| + V.(p.)=0, 
the force equation, 



dv 



V ■ Vw 



-VP 



dt p cp 

the magnetic induction equation. 



dB 

dt 



V X 



V X B 



jxB , {jxB)xB 



c-ypip 



and Ampere's law, 

c 



-V X B. 



(la) 



(lb) 



(Ic) 



(Id) 



Our notation is standard: p is the mass density, v is the velocity, P is the gas pressure, B is the magnetic field, and j is the current density. 
The density, velocity, and pressure all refer to the dominant neutral species. The combination 7pi is the neutral-ion collision frequency, with 



7 



mi 



being the drag coefficient, and jic the electron number density. The quantity (aw) in is the average collisional rate between ions of mass mi 
and neutrals of mass nin', it is equal to 1.69 x 10^^ (T/IO K)^''^ cm'^ s^^ for HCO+-H2 collisions, and is almost identical to this value for 
Na+-H2 and Mg+-H2 colUsions (see McDaniel & Mason 1973). 

The three terms on the right-side of Equation fTct represent induction, the Hall effect, and ambipolar diffusion, respectively. The 
difference between ion and electron velocities gives rise to the Hall effect, whereas t he difference between ion and neutral velociti e s give s 
rise to ambipolar diffusion. Discussions of their relative magnitudes can be found in Balbus & Terquenil (2001) and Sano & Stone (2002). 
For the sake of completeness, however, we repeat here the relative ratio of the ambipolar to Hall terms: 
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where n is the number density of the neutrals and Cs is the isothermal sound speed. Assuming that the final two factors are each about 0.1, 
we see that a neutral density below about 10® cm~"^ brings us safely into the ambipolar diffusion regime (see fig. 1 of Kunz & Balbus 2004). 
Since we are concerned here with the interplay between non-ideal MHD effects and shear, we ignore Ohmic dissipation, for which no shear 
instabilities are present. In writing Equation i fTct . we have implicitly assumed that the thermal-pressure force on the ions and electrons is 
negligible compared to the electromagnetic and coUisional forces. This is an excellent approximation for the systems of interest. In addition, 
the inelastic momentum transfer by the ion and electron fluids due to attachment onto grains and neutralization is negligible compared to the 
momentum transfer due to elastic collisions, and it has been implicitly omitted from the induction equation. We have also ignored coupling 
via ionization and recombination, since these processes are slow compared to elastic processes. 

Consider a flow in a stationary ('lab') Cartesian coordinate system along the y-axis, v = v{x)e.y. Although we take the density and 
the magnetic field to be everywhere uniform, we put no restrictions on the x-dependance of the velocity field and the orientation of the 
magnetic field. Eulerian perturbations to this flow are allowed, denoted by a 5. Keeping only terms linear in 5 and working in the Boussinesq 
approximation, we find 

V ■5v = Q, (3a) 
dSv , ,dSv , , dv(x) „fSP B-SB\ 1 ,„ 

^ ^^^^dSB _ ^ (5 . _ / cB^\ + f }IA21X] V . X SB) + ^VHB , (3c) 



dt dy dx \ ATveric J \ jpi J jpi 

where 

B 

''^ ^ (47rp)i/2 
is the Alfven velocity. 



2.2 Shearing Sheet FormaHsm 

The equations are first transformed from our lab-frame coordinate system to one comoving with the flow, centered at a fiducial location 
(xo, 2/0, zo) moving at velocity v{xo) = vo- We then consider a local neighborhood surrounding this point and Taylor expand the velocity 
field about xq to find 

, , , dv 

v(x) ^ vo + (x — xo) — 

As is weU kno wn, a shearing background precludes simple plane wave solutions to the perturbation equations given in the previous section 
I Goldreich & Lvnden-BelL1965.) . This difficulty is circumvented by adopting shearing coordinates, given by 

x' = X , (4a) 
y'^y-2Axt, (4b) 
z ^z, (4c) 
t' = t, (4d) 
so that 

^ = A _ 2At'^^ , (5a) 
ox ax' oy' 

^ = 7r7' (5b) 

oy ay' 
oz oz' 

where we have defined 2A = (di'/da;):rg01n this frame, the velocity field is v{x') — 2Ax' . The benefit of this coordinate transformation is 
that a spatial dependence exp[i{k'^x' + kyy' + k'^z')] may be assumed for the perturbations, so long as we replace a fixed x wavenumber 
with a shearing one: 



^ The parameter 2A used in this paper is not to be confused with the Oort A constant, which concerns rotating systems. Here, we are primarily interested in 
planar shear flow, and 2A is to be identified as the characteristic frequency associated with the velocity profile. 
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^ k^{t) = fci - 2Ak'yt' . (6) 

No modification of y' and z' variables are needed. Enacting tliis transformation and dropping tire primes for ease of notation, our equations 
become 

kxSvx + kySvy + kzSvz = , (7a) 
+ ikx{ — + — -I- SBx^O, (7b) 



dt \ p 47rp J Aixp 

dSvy , , , fSP B-SB\ fc-B 

+ 2ASvx + iky I h — - I- SBy = 0, (7c) 



dt VP 47rp J Aixp 

dSv, ., fSP B-SB\ k-B 

+ tk^( h — - I- (5Bz = , (7d) 



dt \ p Aixp J Amp 

'^^^^ vju T3\x c(fc-B) ^ c„s , k^VA {kXVA)x N r, /T ^ 

— life • B)dVx (fe X SBjx H SBx - — (fe X VA) - SB = 0, (7e) 

dt ATveric jpi jpi 

- 2^155, - i(k ■ B)5vy - (fe X 5B)y + ^ SBy - ^ ^"^^^ (kxvA)-SB = 0, (7f) 

^_,(fe.B)5„_4^(fex5B). + ^5B.-i^^i^(fex«A)-5B = 0, (7g) 
at 47rene jpi 7Pi 

where the perturbations are now time-dependent Fourier amplitudes and we have suppressed the explicit time-dependent notation in k^ and 
k^ = kx + ky + k'^. Since dkx/dt = —2Aky, Equations fTet - \lg\ together with Equation fTat guarantee the divergence free condition 
d{k-5B)/dt = 0. 

In an attempt to keep the presentation as simple as possible, we first consider Equations (7) in the limit where ambipolar diffusion is the 
dominant non-ideal MHD process (fj3]l, ignoring the Hall effect for the time being. Not only does this aid in our interpretation of the physics, 
but also the two processes generally act in distinct regions of parameter space. We then isolate the Hall effect in 21 neglecting ambipolar 
diffusion. The similarities and differences of the two resulting instabilities are discussed in Sections |5] and [6] Much of the formalism for 
understanding the Hall-shear instability is developed in the following section on ambipolar diffusion, and so it behooves us to encourage any 
readers primarily interested in Hall physics not to bypass the following section. 



3 AMBIPOLAR-DIFFUSION-SHEAR INSTABILITY 



In this Section, we are concerned solely with the interplay of ambipolar diffusion and shear, and we neglect the Hall terms in Equations 
(I7el>-j7g|). Before we begin reducing these equations to a more manageable set, however, it is of interest to note that, while the presence of 
shear causes the y-component of the background magnetic field to grow linearly with time: 

Byit) = By{Q) + 2ABxt, 

where By (0) is the initial y field, the combination k ■ B is constant with time, despite the fact that neither the Eulerian wavenumber 
fe = [kx{t), ky, kz] nor the magnetic field vector B is individually constant: 

fe • B = k'xBx + kyBy{0) + k,B, . 

Unfortunately, the same does not hold for the combination k X B, and so the ambipolar diffusion terms in Equations (l7et-(|7g|( are intrinsi- 
cally time-dependent. This complicates matters. 

We simplify the set of Equations (7) as follows. Using fc • SB = 0, we first obtain 

1 



SB, = 



k. 



kxSBx + kySBy 



Then we may eliminate SB, from Equations fTet and ( iTft to find 
1 



^ + k^Vxx ) SBx + k^rixySBy 



^''-~-^{k.B) 

Here we have introduced the resistivity tensor rj, whose elements are given by 



— + k^riyy j SBy + (^^T^yx " 2A^SBx 



7Pi 



(fc X VA)i (fc X VA)j kj (fe X VA)i (fe X Va) 



7Pi k, 7pi 

where Sij is the usual Rronecker delta function and fe is the unit wavevector. Next, we rearrange Equation fTct : 



(8) 

(9a) 
(9b) 

(10) 
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5P_ B SB _ 1 
p Anp ky(k-B) 



d^SBy . d 



df^ 



+ ^ (^k^riyySBy + k^T^y^SB^ + (fe " V a)^ S By + 2A(^k'^T]xy5By + k'^rixxSBj:^ 



(11) 



Inserting Equations ( |9a| ) and (TTJ into Equation ( I7bb . we obtain 



dt2 



+ ^ ( k'n^^xSBx - -j^k'rtyJB, ) + (fe • i;a)^(5B:. - 2A^ ( k'-q^^ + ^fc^r?^, ) SB.. 



kx 

ky 



d^SB. 



dt^ 



- + 4: [k'^nVySBy - tlk^'ri^ySBy ) + (fc • VA^SBy + 2 A^ k"" T, ^ y 5 B y 



dt 



: krr 



h.2 



(12) 



For economy of notation, we iiave defined k^y = k^ + ky. We need anotlier independent differential equation coupling 3Bx and SBy. 
Multiplying Equation ( iTct by k^ and Equation l l7dt by ky, then subtracting one from the other, we find 



, dSvy „ . , „ , dSvz ., k ■ B ., fc-i?-„ 

kz — z-^ + 2AkzOVx — ky — ; 1-kz— dBy + iky— oBz = . 

dt dt 47rp Anp 

Using Equations fTat . l[8j, and l |9bt in Equation l ll3t leads after some simplification to 



(13) 



d^SBa: 



df^ dt 



+ ^ ( k''rixxSBx 



^k^y^SBA - AA^^ + (fe ■ VAfSB, + 2A (t^—^k%^ - ^k^yA SB, 

+ I (f^VySBy + t^k'v^ySBy^ + (k ■ VA^SBy + 2A^k'r^,ySBy 



kx ky 



(14) 



Again, for economy of notation, we have defined fc^^ = fc^ + fc^. 

For numerical work, it is convenient to isolate the second-order time derivatives. Equations M2\ and ( I14t may be recombined to yield 



(fSB^ 
dt2 



d_ 

dt 



{k^rixxSBx + k'^T]xySBy^ + AA 



d^SBy d 



— [k'^'qyySBy + fc^T^y^^B^^ + 4:A 



kx ky dS X 

k2 dt 

ky dSBx 



(fc • VAfSBx + 4:A^ {k'r^xxSBx + k^n^ySBy^ , 
[k^r)xxSBx + k'^r)xy5B^ 



(fc - VAfSBy - 2A y 



Equations l |15l l and 1 II6I 1 are the two coupled differential equations in SB, and SBy that form the cornerstone of the analysis. 



(15) 
(16) 



3.1 Qualitative Behavior 

In a system where ambipolar diffusion and shear are absent, it is straightforward to show from Equations il5\ and il6l that the perturbed 
magnetic field lines simply follow fluid elements, resulting in Doppler-shifted Alfven waves propagating along the background magnetic 
field. The introduction of shear into the picture has two effects. The first effect is that vorticity is generated in the flow. The accompanying 
centrifugal force (associated with the resulting eddy) pushes on the shear interface, resulting in the growth of any deformation in the interface 
with time, provided that the vorticity is out of phase with the surface deformation. This is the essence of the Kelvin-Helmholtz instability. 
The second effect is that any x-displacement in the magnetic field becomes sheared out into an j/-displacement. Thus, the perturbed magnetic 
field is effectively rotated until it becomes aligned with the shear interface. 

In the presence of shear, any physical mechanism that conspires to rotate SBy back into SB, completes a feedback loop and results in 
growth. Ambipolar diffusion does just that. Since ambipolar diffusion only affects those currents flowing perpendicular to the background 
magnetic field, it tends to align magnetic field perturbations perpendicular to the background magnetic field. This manifests itself as an 
effective rotation of SBy into SB, (albeit with a decrease in |5iJ|). In this case, ambipolar diffusion and shear conspire to stretch any 
perturbation in the magnetic field, resulting in an exponentially-growing instability. In a gas where either (1) the magnetic field is so strong 
that its tension effectively resists being stretched by the shear [i.e., (fc • va)^ 3> 4:A^], or (2) the bulk neutral fluid is so poorly coupled to 
the magnetic field so that its velocity does not grow with it [i.e., {-ypi) <S 2A], this instability does not operate efficiently. We will show, 
however, that there still remains a great deal of unstable parameter space with which to work.^ 

This route to instability has been seen before in the ambipolar-diffusion-modified MRI IKuuz & Balbusll2004 ; Deschll2004 ), where the 
role of shear is played by the differential rotation of an accretion disc [v{x) = xQ{x), where fl is the orbital frequency]. Here, however, 
we see that this destabilizing behaviour is part of a more general process, and in no way depends on rotational kinematics. The finding that 
ambipolar diffusion renders an accretion disc unstable (albeit weakly) for both inwardly- and outwardly-decreasing angular velocity profiles 
acquires a natural explanation in the more general context of this paper. Ambipolar diffusion can destabilize any shear flow profile, in very 
much the same way that the sign of dv{x)/dx does not determine the outcome of the Kelvin-Helmholtz instability. In fact, we will show that 
the criterion for this ambipolar-diffusion shear instability is independent of the magnetic field strength, and is reliant only upon the ratio of 
the ion-neutral collision frequency 7pi to the frequency implied by the shear of the flow, 2A, and the geometry of the magnetic field. 



3.2 The Case ky = B, = 0: & Time-Independent Zero-Order State 

The qualitative behaviour discussed above is most easily seen in the simple case of ky = B, = 0. In this situation, the zero-order state is 
time-independent, and we may seek solutions to Equations dlSt and l ll6t with time dependence exp((jt). The resulting dispersion relation is 



6 M.W. Kunz 



k^vl + (k-VA) 



7Pi 



(T + Co = , 



where 



Ca = 2{k-VA) + 



k^v\{k ■ va) 



2A 



k^By{k-B) 



Co = (k- va) 



(k ■ va) 



2A 



k^By{k-B) 



Before proceeding to obtain an instability criterion, let us note that Equation ( |17| l may be written in the more compact form 



7 2 2 

fc '"A 

7Pi 



cr + (fe ■ va) 



+ 



{k-VA? 

7Pi 



CT + (fe • IJa) 



2yl 



k,By{k-B) 

inypip 



a +(fe-i»A) 



(17) 



(18a) 



(18b) 



(19) 



Consider the limit of vanishing shear. In this case, the right-hand side goes to zero and the two brackets on the left-hand side become 
decoupled from one another. The two solutions obtained by setting the left bracket to zero correspond to forward- and backward-propagating 
Alfven waves with 55 || B, which are damped at a rate k'^v\/2')pi. The other two solutions obtained by setting the right bracket to zero 
correspond to forward- and backward-propagating Alfven waves with SB || (fc X B), which are damped at a rate (fc ■ va)'^ /2^pi. One 
consequence of the difference in these damping rates is an effective rotation of SB ■ B into SB ■ (fe X B). As discussed in i]3.1l this 
difference is at the heart of the ambipolar-diffusion-shear instability. 

A sufficient condition for unstable solutions to exist in Equation l ll7t is Co < 0, or 



(fe ■ va) 



2A 



k^By{k-B) 



< 0. 



(20) 



This is identical to the result in §4.2 of Kunz & Balbu j i 2004 ) in th e limit of vanishing rotat ion frequency Q, (so that the epicyclic frequency 



K. 0), and is similar to the instability criterion for the ideal MRI JBalbus & HawlevI 19911) : 

k1 dn^ 



Notice, however, that the stabilizing effects of ambipolar diffusion seen in the second term of equation (35) of Kunz & BalbusI ( 2004 ) [which 
is oc ii^{k ■ V a)"^ / {') pij^} are absent here. This term represents epicyclic oscillations in the neutral fluid, an inherently stabilizing motion 
in Keplerian discs, being communicated to the (potentially) unstable ions through collisional coupling. In the limit 7Pi — > (i.e., infinite 
neutral-ion collision time-scale), it is clear that this term dominates and the bulk neutral fluid oscillates at the epicyclic frequency, unaffected 
by the presence of the magnetic field and the ions that are tied to it. In other words, if the neutrals can respond to magnetic forces on an 
epicyclic time-scale, the MRI will act on both fluids and will be effective in transporting angular momentum. Since it is rotation that gives 
birth to these stabilizing epicy clic oscillations, th is term is absent in the dispersion relation given here. Evidentally, the situation investigated 
in both Kunz & BalbusI 1 20041) and DeschI 1 2004 ) was the superposition of two different, but related, instabilities acting in tandem: the MRI, 
whose magnetic 'tether' between fluid elements, essential to the transport of angular momentum, is undermined by the imperfect coupling 
between the ions and neutrals (leading to decreased growth rates); and an ambipolar-diffusion-shear instability, which results from the 
combination of shear and the anisotropic damping of ambipolar diffusion. 

Further insight into the interpretation of Equation l l20t is afforded by rewriting it in the form 



2A ^ k^B^ 
7pi k^By 



(21) 



Note that the strength of the magnetic field is not at all relevant here; only the ratio of the shearing time-scale to the neutral-ion collision 
time-scale and the geometry of the magnetic field come into play. Furthermore, the freedom in choosing the sign of kz/kx guarantees that 
any non-constant velocity profile can be destabilized, regardless of the sign of its derivative. Physically, this equation states that the time for a 
neutral particle to collide with an ion must be longer (by at least the factor given on the right-hand side) than the time it takes for a magnetic 
perturbation to grow by shear. If this condition is not met, the neutral fluid is well-coupled to the magnetic field, and we are left with simple 
linear-in-time growth due to shearing of the magnetic field perturbation. 
Defining the dimensionless parameters, 

^^(k2VAf_ Am=^, (22) 



4^2 



M 

7Pi 



the dispersion relation may be written in dimensionless form and growth rates may be determined numerically. In Fig. [T] we give three- 
dimensional plots of growth rate in the X-\kx/kz\ plane (for Am — 1) and in the Am-\kx /kz \ plane (for X — 1). The signs of kx/kz and 
2 A are chosen such that instability is possible, and By/Bz is taken to be unity. Increasing By/Bz does not significantly affect the growth 
rates, but rather opens the available unstable space to smaller values of Am. Note that there is less unstable parameter space as one goes to 
small Am (the fluid becomes well-coupled to the magnetic field). The boundary separating stability from instability is given by Equation 
( I21b . The maximum growth rate (~ 0.1 \2A\\ see Equation [33]below) is shown in Fig.[T}; for the parameter space spanned by X and Am. 
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Figure 1. Counter-clockwise from top: ambipolar-diffusion-shear instability growth rates in (a) tlie X-\kx/k^\ plane with Am = 1 and (b) the Aia-\kx/kz\ 
plane with X = 1; (c) maximum growth rates in the Am-X plane. All figures have By/B^ = 1. Only regions of instability ai'e shown, with the height being 
proportional to the growth rate. The maximum growth rate is ~ 0.1 |2A| (see Equation 133) . Note that the boundary separating stability from instability is 
given by Equation j2U . 



One final comment is worth mentioning concerning Equation l l21t . In a time-dependent zero-order state {ky, Bx 7^ 0), the denominator 
of the right-hand side grows quadratic in time and the instability criterion will quickly become trivial to satisfy. It is therefore of interest to 
rigorously test whether this situation is realizable by performing an analysis of this more general case. 



3.3 The General Case 

We now investigate the general case of mixed wavenumber and field geometry via two approaches. First, we employ a WKB technique. While 
this approach is strictly applicable only under special circumstances, it will aid in the interpretation of results from our second approach: a 
direct numerical solution of Equations l ll5t and l ll6t . 



3.3.1 WKB Analysis 

Without loss of generality, one can express each perturbation variable in a WKB form: 

5{t) = 5'"\t) exp fi f oj{t)dt 

n=0 ^ 

where the WKB phase has been expressed in terms of an integral over a slowly varying frequency. This is valid as long as <^ (j'"' 

and dXxiLo/dt <^ lo (the adiabatic approximation). Both of these conditions are guaranteed if ky/k^ <g 1, and we may immediately identify 
the WKB parameter k^/ky. Physically, the WKB parameter represents the ratio of the time for k to change significantly to the shearing 
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time-scale[f| Expressing SBx and 5By in tlie form i23t . taking the limit ky/kz 0, and defining a = iuj, the lowest-order terms from 
Equations il5l and il6l result in the following dispersion relatiorlfl 



a'^ + k^ tr(77)a^ + C2cr^ + k^ tv{r]) (k ■ va^cj + Co = , 
where 

Ca = 2(fe • VA.f + k^ det(j7) + 2Ak'^r)^y , 
Co = (fe • va)^ ^{k ■ VA f + 2Ak'^r]a:y^ . 



(24) 

(25a) 
(25b) 



Here, tr(j7) and det(?7) denote the trace and determinant, respectively, of the resistivity tensor rj. This is similar to the dispersion relation 
([nil in El and may be written in a form similar to that of Equation l ll9t : 



2 I 



where 



+ (k-VA) 



2 I 



fc^tr(77) 



+ (k-VA) 



1/2 



~2Ak'^rixy l^o-^ + (fe ■ ua)^] , 



(26a) 



(26b) 



Evidently, the effect of shear is to couple different polarizations of Alfven waves, which are damped at different rates. It is also of interest to 
calculate the associated eigenvectors. Defining 



5-^ = CT^ + k'hr{ri)a + k"^ det{ri) + (fe ■ VAf + 2Ak'^r)^y , 

the (non-trivial) eigenvector components (in the limit ky/k^ <C 1) can be expressed as 

5vx_ 

^ aa^ + k'^r/xxik-VA)'^ 



^~2Aa^ + (fe ■ VA)^k'^Jiyx^ 



^—^^^[k\^a + 2Ak%.], 
i{k ■ B) Svx 



\a'^ + k^rjyyG + (fc • wa)^ + 2Ak^rixy 



o-CT^ + k'^l]xx{k ■ vaY 
As in i]3.2| a sufficient condition for instability is Co < 0, or 

{k-VAf + 2Ak^'qxy < 0. 



(27) 

(28a) 
(28b) 
(28c) 

(29) 



This criterion depends on both kxit) and By{t), through -qxy As a result, the destabilizing term will grow in amplitude, opening up more 
and more unstable parameter space as time progresses and increasing the growth rate. 

If we view the dispersion relation ([24} as an equation in k^ and (fe ■ ua)^, the maximum growth rate can be calculated for a given rj. At 
the maximum growth rate a — (Jmax, partial differentiation of Equation ( I24t with respect to k^ and (fe • va)^ gives the two equations 



fc^tr(77)crmax + [2/0"* det(r7) + 2Ak^rixy'^ (Jmax + fc^tr(T7) (fe • i;A)^crmax + 2Ak^r)xy (fe • wa) 



0, 



2a' 



fc^tr(T7)crmax + 2(fe • VAf + 2Ak''rjx 



Eliminating (fe • vaY between these two leads after regrouping to a surprisingly simple result, 

[4det(j7) - tr2(j7)jo-Lx - 4Ar7,ytr(T7)crn,ax - 4^%^^ = . 

There are two solutions to this equation, corresponding to the extrema of the dispersion relation. Only one of these is a physically meaningful 
solution satisfying the dispersion relation it is given by 



(30) 
(31) 

(32) 



1^1 



"Tjxy 



deti/2(r,) + itr(T7) 



(33) 



^ There is actually an additional condition we must impose. The resistivity tensor 77 is dependent upon fc X da, and thus has an intrinsic time-dependence. 
We therefore require that fc X va is a slowly-vaiying function of time, or equivalently, that By/Bx 2> 1- Such a lai'ge By trivially satisfies the instability 
criterion. 

^ Retaining terms of order ky/kz contributes an extra term, similar to the final term in equation (2.25) of lSalbus & HawlevI jl992lj) . that does not affect the 
essentially oscillatory or exponential behaviour of the solution; it is —iAkxkyk"^ [cro-^ + k'^rjxxik • VAf], with given by Equation Its importance 
in a WKB treatment is as an amplitude modifier on longer time scales. In a stable system, it represents the competition between shear, which is trying to stretch 
the perturbations to result in linear-in-time growth, and ambipolar diffusion, which is trying to dampen this growth. 
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Figure 2. Evolution of the perturbed magnetic fields 53^ (solid line) and 5By (dashed line) for (a) (X, Am) = (1.0, 1.0) , (b) (X, Am) = (1.0, 0.1), (c) 
{X, Am) = (0.1, 1.0), and (d) (X, Am) = (0.1, 0.1). The sign of 2Ak'^/kz is chosen to be negative so that the flow is stable. The initial amplitudes are 
5Bx = 0.3 and SBy = 5Bz = 0, and ky/k^ = Bx/B^ = By{0)/Bz = 1. The initial value oi t = k'^/kz is determined by k • SB = 0. 



This is a remarkable result. ^B albus & Hawlev 1992a ) conjectured that the maximum growth rate of any instability feeding off the differential 
rotation in a disc is given by the local Oort A value, a a = {l/2)\dQ/d In _Rj, no matter the cause for instability. The reason is rooted in the 
dynamics of the differential rotation itself. In this paper, we are concerned with planar shear flows, and we arrive at a similar result, with the 
shear playing the role of the differential rotation. It is notable that Equation ( 133b is independent of the degree of ionization, depending only 
upon the geometry of the background magnetic field and the shearing rate \2A\. 

The next order in a WKB expansion yields the time dependence of the slowly- varying amplitude. Provided we are able to compute this 
amplitude, the eigenvectors may be used to compute the wave energy and compare with the numerical results given in the next section. This 
has been done, e.g., by Johnson (2007) for the case of nonaxisymmetric shearing waves in differently-rotating disks. In a non-dissipative, 
continuous system, the amplitude may be computed from conservation of wave action. In the presence of ambipolar diffusion, however, 
wave action is not conserved, and retrieving the amplitude in this fashion is prohibitive. Instead, one must take the algebraically tedious 
approac h of directly ca lculating the higher-order expansion for the modes. We have done this, finding an equation similar to that of equation 
(A5) of Johnsorj (2007). Unfortunately, once the eigenvectors ( I28ab - J28ct are substituted in, the result cannot be easily integrated and the 
slowly- varying amplitude cannot be calculated. Numerical solutions seem to be the most profitable approach. 



3.3.2 Numerical Solution 



Here w e undertake a direct numerical solution of Equations dlSb and ( II6b . Following lines similar to those developed in lGoldreich & Lvnden-Beli 



( 1963), we introduce a new independent time variable 

_ kx{t) _ kx r, A^y , 



(34) 



kz kz kz 

so that {k/kzY ~ t"^ + (ky/kz)^ + 1. Equations dlSb and ( II6b may then be written in dimensionless form and numerically integrated. 
All that remains is to specify initial values for SB^, SBy, 5Bz, and By/Bz. (The initial value of r is determined from k-SB — 0.) 
Unfortunately, we must also specify values for ky/kz and B^/Bz. However, we have undertaken a parameter study to see if varying 
these significantly influences the results, and have found that the qualitative behaviour is not affected. The results presented here have 
5Bx{0) = 0.3, 5By{Q) = 5Bz{<d) = 0, and ky/kz = B^/Bz = By{0)/Bz = 1. 

In Fig.m we give the evolution of SB^ (solid line) and SBy (dashed line) for (a) {X, Am) = (1.0, 1.0) , (b) {X, Am) = (1.0, 0.1), 
(c) {X, Am) = (0.1, 1.0), and (d) {X, Am) = (0.1, 0.1). Here, the sign of 2Ak'x/kz is chosen to be negative so that the instability does 
not operate. The modes are damped due to ambipolar diffusion at a rate proportional to 2A/'ypi (i.e., a longer neutral-ion collision time-scale 
results in faster ambipolar diffusion). Decreasing X results in a smaller damping rate (i.e., a longer Alfven crossing time-scale allows more 
time for an ion to communicate the presence of the magnetic field to a neutral via collisions). 

In Fig. [3] we give the evolution of the corresponding unstable modes for the same parameters as in Fig.|2] The evolution of 5 By in this 
case is similar to 5Bx, and therefore is not given. Note that In \5Bx\ is plotted, rather than SBx, so that the growth may be more clearly 
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Figure 3. Same as in Fig.|2] except that the sign of 2Ak'^/kz is chosen to be positive so that the flow is unstable. Note that the ordinate value is In 
rather than SBx as in Fig.|2] The evolution of SBy is similar and therefore not shown. There are, in general, two stages of evolution. Before the instability 
criterion is satisfied (i.e., t = kx{t)/kz is too small) we have wave damping due to ambipolar diffusion. Once r = kx{t)/kz reaches the value where 
Equation (29) is marginally satisfied, the growth phase begins. 



seen. There are, in general, two stages of evolution. Before the instability criterion is satisfied (i.e., t = (t) /fez is too small) we have wave 
damping due to ambipolar diffusion (i.e.. In \5Bx\ < 0). (The cusps on the plots correspond to the wave amplitude passing through zero and 
changing sign. For example, in Fig.[3]3, the wave executes one full cycle before becoming unstable.) Once r — kx{t)/k2 reaches the value 
where Equation l |29l ) is marginally satisfied, the growth phase begins (i.e., In \5Bx \ > 0). 

Similar to investigations of other instabilities in the presence of a shearing time-dependent background (e.g., Balbus & Hawley 1992b, 
Balbus & Terquem 2001), we find that the evolution unfolds as a series of time-independent problems, whose analysis is aided by the 
existence of WKB solutions. Here, however, the analysis is complicated somewhat by the explicit time-dependence of the unstable terms. 
To keep things simple at first, consider the case Bx = 0, so that the time-evolution of the destabilizing term is dependent only upon the 
wavenumber kx{t){= k'^ — 2Akyt). In Fig.|4^ we show the evolution of an initially leading wavenumber through the parameter space 
defined by Am and \kx/kz\ for ky/kz — 1 and By{T) — B y(0) = 0. The wavevecto r evolution traces (and retraces) a horizontal line in 
the plane of Fig.|4^. In contrast to the nonaxisymmetric MRI Balbus & Hawlevlll992bh , the wavevector evolves towards stability, not away 



from it, and in principle spends only a finite amount of time in the stable region, that is, when (fc ■ va)^ > —2Ak^rjxy is satisfied. In Fig. 
|4}5, we relax our assumption of a vanishing Bx and take Bx/ B^ = 1, so that By is time-dependent. Recall that kx/k^ is our time parameter, 
so that By increases with fc^. One modification to the picture presented in Fig.|4^ is an extension of the unstable region at large \kx/kz\. This 
is due to the time-dependence of By (and thereby, rjxy)- Another important consequence is that, in principle, the growth rates will increase 
in time due to the fact that a strongly trailing (or leading) wavevector is trivially unstable. (These modes are well-localized in a WKB sense.) 
If the exponential growth phase goes unchecked by non-linear effects, it is conceivable that the growth rate will quickly surpass the shearing 
rate, or for that matter, any other dynamical rate in a system of interest. 

One final issue remains: where does the energy come from to sustain this growth? The last term in Equation l |29t introduces a novel 
form of coupling: in dyadic notation, it is 

k^ri : V I) . 

This suggests that ambipolar diffusion influences the free energy path between the velocity shear and the perturbations. In an ideal MHD 
shear flow, where ambipolar diffusion is absent, the link between the fluctuations and the free energy source is provided by vortex stretching. 
Magnetic fields modify this picture, of course, providing what amounts to an effective surface tension by resisting the shear. The effect of 
ambipolar diffusion may be most easily seen by using Equations l llSt and l ll6t to derive an equation for 

SJz = kxSBy - kySBx , (35) 
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Figure 4. (a) Regions of stability and instability in the Am-\kx/kz\ plane for Bx = 0. Initially, \kx/kz\ is large, and the point defining the wavevector 
moves to the left through the unstable region on a constant Am line. If (fc • da)^ > —2Ak^r]xy, then a finite portion of time is spent in the stable region. 
After attaining its minimum (i.e., k^/k^ = 0), the wavevector point retraces its path to the right, reentering the unstable region, (b) Same as in (a), but for 
Bx/Bz = 1. See text for discussion. 




Figure 5. Evolution of (a) SJz for k'^rj : V i! > (stable), and (b) in \5Jz | for k'^rj : V i! < (unstable). Here, {X, Am) = (1.0, 1.0); initial values are 
as in Figures [2]and[3] 



which is proportional to the ^-component of the perturbed current density: 



J2 



|5 + (fc.«Af 



AAky dSBy 

U J z 



^riyySBy + k^rjyxSBx^ + ^k^ T^x X S B x + rjxy 5 B y^ — 2A{k^rixx5Bx + k^rixySB^ . (36) 



dV 

In the absence of ambipolar diffusion, the right-hand side vanishes, and we are left with a wave equation for SJ^ with the additional effect 
of shear acting as an amplitude modifier. The right-hand side is best approached piecemeal. The first two terms represent the damping of the 
current due to ambipolar diffusion. It is clear that different components of the perturbed current are damped differently. If ambipolar diffusion 
acted isotropically, with a resistivity rj = rjxx ~ rjyy, these terms would read 

f\jX U,L t\jx 

which shows the effect of resistivity on the propagation of the wave, as well as the attempt to resist the stretching of the perturbed magnetic 
field by shear. The final term in Equation l l36t is the key to instability. It becomes a source term when k^r\ : V < 0, extracting energy from 
the background velocity shear. In Fig.|5] we show the evolution of &Jz for (a) yp'r] : V w > (stable), and (b) iP'r] : V i; < (unstable). For 
the unstable case, the z-component of the vorticity (= kxSvy — kySvx) follows a similar evolution to the current. 



4 HALL-SHEAR INSTABILITY 

We now consider the case where the Hall effect is the dominant non-ideal MHD process, corresponding to a neutral number density much in 
excess of 10^^ cm~^. The analysis follows the same course as in §[3] provided we redefine the resistivity tensor to be 

^ ckz{k-B) ( kxky/kl klz/kl 

GTIq \ ~ ' kxz I kj z ' ky I kj 2, 
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Realizing tliat A:^tr(?7) = and A:* det(j7) — (?k^ [k ■ B)^ / (inenc)'^, the WKB dispersion relation (see i]3.3.U may immediately be written 
down: 

cfc^(fc-B) 



2 . ck(k ■ B) ,, 
+ I ^ + (fe ■ va) 



ck{k-B) , 2 

r— h (k-VA) 

4-Kene 



-2A- 



■ p + (fe ■ VAf \ ■ (38) 



This is identical to equation (81) in lBalbus & Tergueml d2001h in the limit of vanishing rotation frequency Q, (and negligible Ohmic dissipa- 
tion). Note that it is not necessary to assume By/ <C 1, as in ij3.3.1l since the only time-dependence is due to the wavenumber k (recall 
that the combination fe ■ i? is a constant). 

Before we proceed any further in analyzing Equation l l38b . it pays to examine the simple case of vanishing shear. The dispersion relation 
then becomes 

2 . ck{k-B) n2 n ^1C^^ 

a ±1(7 — ^ - + (k-VA) =0. (39) 

This is precisely the dispersion relation for a uniformly magnetized plasma with the displacement ciu'rent and electron inertia ignored (see, 
e.g., Krall & Trivelpiece 1973); its positive frequency solutions are 



.2 , c^k^jk-By ^'^' 

(k-VA) H ^ 



ckjk-B'' I" ^2, an. D^2 



g47^2g2j^2 



(40) 



For small wavenumbers (low frequencies), both of the above solutions reduce to Alfven waves. At larger wavenumbers, however, right- 
handed waves (plus sign) go over to the high-frequency whistler wave branch, whereas large k left-handed waves (minus sign) are cut off at 
a frequency 

eB /nc\ _ /nc\ 

t^cutoff = = I^CM , (41) 

fic \ n / \ n / 

where /j, is the mean mass per neutral particle. These modes arise due to the drift of the field lines with respect to the ion fluid as an Alfven 
wave propagates through the lighter electron fluid, and are related to the so-called R- and L-waves of plasma physics. 

Restoring the velocity shear, it is straightforward to show that instability proceeds through the point a = 0. The instability criterion is 
therefore 

(fc.^A)^ + 2A "^-(^"^^ <0. (42) 

47reno 

This may be written in a more physically transparent form: 

2A ^ _k^ /n,^ ^ ^^3^ 



(jJcp kzB \ n 

In other words, the time required for an ion to execute one orbital gyration around a magnetic field line must be longer (by at least the 
factor given on the right-hand side) than the time it takes for a magnetic perturbation to grow by shear. If this condition is not met, the ions 
are well-coupled to the electrons (and thereby the magnetic field), and we are left with simple linear-in-time growth due to shearing of the 
magnetic field perturbation. The freedom in choosing the sign of fcz(fe ■ B) guarantees that any sign of shear can be destabilized, similar to 
the result in 3. 2 1 involving ambipolar diffusion. 

The maximum growth rate of this instability may be calculated directly from Equation ( |33l l: 

O-max = 1^1 , (44) 

which occurs when k = k-^Bz. This is in agreement with the findings of iBalbus & Tergueml 1 2001 ). despite our neglect of rotation. Defining 
the dimensionless parameters, 

v-_ {k-VA? , TT - 1 ck-.{k-B) 

X = — and Ha = - — , (45) 

4yl2 \2A\ A-Kcn^ 

the dispersion relation may be written in dimensionless form and growth rates may be determined numerically. In Fig. |6] we give three- 
dimensional plots of growth rate in the X-\k/kz\ plane (for |Ha| = 2) and in the |Ha|-|fc/fc2| plane (for X — 1). The signs of kz{k ■ B) 
and 2A are chosen such that instability is possible, and ky is set to zero. Note that there is less unstable parameter space as one goes to small 
|Haj (the ion-neutral fluid becomes well-coupled to the magnetic field). The maximum growth rate is shown in Fig. [6}; for the parameter 
space spanned by X and jHa|. It is important to notice that if a strongly-leading (or -trailing) wavevector (|fc/fcz| ^ 1) begins its life in the 
unstable regime, it will always be unstable. By contrast, a stable wavevector will remain stable. 

We had argued earlier (see i]3.1b that, in the presence of shear, any physical mechanism that conspires to rotate SBy back into SBj; 
completes a feedback loop and results in growth. HEMFs accomplish this task not by preferential current damping, as in the case of the 
ambipolar-diffusion-shear instability, but rather by current generation. Consider the case of a purely vertical {z) wavenumber and magnetic 
field. The effect of shear on the background magnetic field and the coordinate system vanishes, and we may write down the linearized 
induction Equations f/et and ( fTft in the form 

d5B^ ckz{k-B) 

— 7^ 1 7^ SBy = i(k ■ B)dVx , (46a) 

at Anenc 
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Figure 6. Counter-clockwise from top: hall-shear instability growth rates in (a) the X-\k/k2 1 plane with |Ha| = 2 and (b) the |Ha|-|A;/fez | plane with X = 1; 
(c) maximum growth in the X-|Ha| plane. Only regions of instability are shown, with the height being proportional to the growth rate. The maximum growth 
rate is \A\. 



dSBy 
dt 



ck,{k-B) 



+ 2A 



5B^ = i{k-B)5vy 



(46b) 



We assume that instability is possible, i.e., that the Hall and shear terms in Equation l l46bt have opposite signs. It is clear from Equation 
l |46b| ) that shear uses SB^ to generate a SBy. The Hall terms, on the other hand, generate SB^ at the expense of SBy. This effect is actually 
present in the absence of shear, and arises because the y-component of the perturbed electron velocity differs from the ion-neutral velocity 
(recall that ambipolar diffusion is ignored here) by a term involving SJy{= k^SB^ ~ k^SB^). The induced magnetic field is sheared further, 
leading to runaway. This behaviour can be seen differently by viewing the Hall terms as 'Coriolis' terms in the magnetic field equations. 
These give birth to magnetic 'epicycles', i.e., circularly-polarized electromagnetic waves. When the sign of the rotation imparted by the shear 
is opposite to that of the handedness of these waves, a struggle ensues over control of the direction of the perturbed magnetic field vector. 
For every bit of field line stretching in the y direction due to shear. Hall forces rotate this increased field back into the x direction, only to be 
stretched further by shear. The compromise of this struggle is an exponentially- growing instability. 



This route to instability has been seen before in the Hall-modified MRI jWardle|[l999l : Isalbus & Terquen]l200lh . where the role of 
shear is played by the differential rotation of an accretion disc [v{x) — x^l{x), where Q is the orbital frequency]. In that case, however, 
the behaviour of the Hall terms is complicated by rotation, which influences the magnetic epicycles by introducing a sense of helicity (see 
section 3 of Balbus & Terquem 2001 for a discussion). The finding that HEMFs render an accretion disc unstable for both inwardly- and 
outwardly-decreasing angular velocity profiles does not depend on rotational kinematics. Rather, it is a simple result of the influence of shear 
on the propagation of circularly-polarized electromagnetic waves. 

At this point a few questions naturally emerge. First, why does the Hall-shear instability require a vertical (z) wavenumber in order to 
operate, whereas the ambipolar-diffusion-shear instability does not? Simply put, the interaction between shear and the Hall effect is strongest 
when the motions implied by the shear lie in the same plane as the magnetic "epicycles" induced by the Hall effect. In other words, instability 
is maximized when the vorticity F = V X the wavevector k, and the magnetic field B all share a mutual axis. Translated mathem atically, 
this im plies that the relevant Hall term is (fc • B){k ■ F); it must be negative for destabilization. This is distinct from the findings of lWardle 
(199i) and lSalbus & Tergueml J200lh . who found that the relevant coupling parameter is, respectively, (f2 • B) or (fc • B){k ■ fl) where 
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f2 is the rotation vector. This corroborates our finding that the Hall instability is not a result of differential rotation per se, but rather of 
any source of shear. Our coupling parameter matches theirs when the vorticity shares the same axis as the rotation, a situation typical for 
disc systems. Second, why is the growth rate for the Hall-shear instability typically so much greater than for the ambipolar-diffusion-shear 
instability? Both processes involve extraction of free energy from the background shear flow. The reason lies in the fact that the Hall effect 
employs a conservative process (cyclotron gyrations) rather than a dissipative process (ambipolar diffusion) to rotate 5 By back into 6Bx- 
Not surprisingly, the difference in the growth rates is related to the rate at which the gas is heated due to ion-neutral friction. 



5 DISCUSSION: SHEARS, ROTATIONS, AND PROJECTIONS 



Despite the impression one may get from the abundance of mathematical manipulations in the preceding sections, the instabilities themselves 
are actually quite simple. They are the result of combinations of shears, rotations (in the case of the Hall effect), and projections (in the c ase 
of amb ipolar diffusion). Much in the way that the MRI may be understood by two orbiting masses connected by a spring (Balbus & Hawleyl 
1992ar) . there exists an equally intuitive toy model that captures the essence of the shear instabilities examined in this paper. Through simple 
matrix multiplication, several main results can be recovered without recourse to the lengthy mathematical manipulations undertaken in the 
preceding sections. All we require is some linear algebra. 

Consider the following abstract problem, in which a position vector |ro) = |2:o, yo) (in Dirac notation), having the coordinates [xq, yo] 
in a two-dimensional x-y Cartesian coordinate system, is subjected to various shear, rotation, and projection operators given by 





1 




s = 






e 1 





cos U — sm t 
sin 6 cos f 



and 



sm <p 
- sin (j> cos (h cos 



sm cp cos ( 

2 1 



(47) 



respectively. We denote by the column vector \r„) — \x„,yn} the position vector |r) after n transformations have been applied to |ro). To 
be precise, the combination S|r) results in a shearing of \r) along the y-axis by a distance ex; R|r) results in a counter-clockwise rotation of 
|r) through an angle 9; and, P|r) projects |r) onto the unit vector | — sin (p, cos (p). These will be put in a physical context below. 

We first turn our attention to the shear and projection operators. Applying the combination SP (a projection followed by shear) to the 
initial state vector |ro) advances it to |ri):0 



|ri) = SP|ro) = (xo sin (j> - yo cos < 



sm (p 
e sin (h ~ cos c 



(48) 



A graphical depiction of this process is given in Fig.|7] While jro) is clearly not an eigenvector of the operator SP, it is straightforward to 
show that |ri) is an eigenvector, with eigenvalue (1 — e sin (p cos (p). It is then trivial to write down the result for general n 1: 



\rn) = (1 - £sin<^cos<^) \rn-i) = (1 



e sm cos < 



(49) 



Put differently, after this transformation has been performed just once, all subsequent applications simply stretch (or contract) the vector 
along |ri) by a factor equal to the associated eigenvalue. Note that in order for any evolution to occur, (f> ^ mn/2 with m an integer. For 
growth, we require 



£ sin (h cos S < . 



(50) 



We now place these results in the context of the ambipolar-diffusion-shear instability. With e = 2AAt and \r) — SB, our shear operator 
corresponds physically to the production of SBy due to the shearing of SBx in a time At. When (p is taken to be the angle between a 
background magnetic field and the x-axis, the projection operator embodies ambipolar diffusion acting on the vector |r) = SB. Their 
combination leads precisely to the behaviour described in i]3.1| In fact, it embodies the exact solution for the case k — kzB^ and B — 
B cos (pBx + B sin (p^y, as can be readily verified by taking the limit At ^ in Equation ([49} to arrive at the differential equation: 

dr 
dt 

If A sm{2(p) < 0, we have exponential growth. The growth rate is A sm{2(p), which has its maximum at precisely the Oort A value. That the 
maximum growth rate occurs at = 7r/4 is in accord with the well-k nown "sin 21" law (e.g., Mihalas & Binney 1981), a result of particular 
significance; we refer the reader to §2.4 of lBalbus & Hawlevi ( Il992ah for a full discussion of this topic. 

Next we turn our attention to the shear and rotation operators. Applying the combination SR (a rotation followed by shear) to the initial 
state vector [ro) advances it to |ri): 



-Asm{2(p) r . 



(51) 



ri) = SR|ro) — {xo sin 9 ~ yo cos t 



+ {xo sin 9 + yo cos ^ 



(52) 



^ This requires \ro) ^ | cos cf>, sin cf>), which con'esponds to a zero eigenvalue of the operator SP. In this case, the projection gives the null vector, and there 
is nothing left to shear For this particular initial state vector the combination PS works fine. 
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Figure 7. Evolution of the initial state vector |ro) (solid arrow) under the transformation SP. (a) A projection is applied to |ro) that retains only its component 
lying along the unit normal | — sin 0, cos </>). The result is P|ro) (grey arrow), (b) This vector is then sheared along the y-axis into |ri) (open an'ow). This 
process represents the ambipolar-diffusion-shear instability (see text for details). Figure is drawn to scale for comparison with Fig.[8] 



A graphical depiction of this process is given in Fig. [8] For ti J? 2, it is possible to show that a recursion relation exists between successive 
state vectors: 



|r„) - (2cos6' - esin6')|r„_i) + |r„_2> = 
This equation has the general solution: 





- A'l\ 






-A_j 





where 
A+ — I cos 9 — — sin ( 



A+ - A- 



cos U sm ( 

2 



1 1/2 



(53) 



(54) 



(55) 



are the characteristic roots of Equation ( 153b . With 9 — uAt and \r) — SB, this transformation corresponds to the physical behaviour 
described in ij4]concerning the evolution circularly-polarized electromagnetic waves (with frequency uj) in the presence of shear; namely, in 
a time At, the perturbed magnetic field vector is rotated by R through an angle uiAt and sheared by S along the j/-axis. Taking the At ^ 
limit of Equation l l53t gives the differential equation 



dt2 



= ~uj{uj + 2A); 



(56) 



(57) 



whose solutions are exponentially growing if 
^<-l. 

UJ 

The similarity between this instability criterion and Equation ( 143 1 is striking. When |j4/a;| equates to unity, the growth rate attains its 
maximum value: precisely the Oort A constant. This simple model captures all the salient features of the Hall-shear instability. 



6 SUMMARY 

In this paper, we have investigated the stability of wea kly-ionized, magnetized planar sh ear flows to linear disturbances. Employing a local 
approach similar to the shearing-sheet approximation of^ Goldreich & Lvnden- Bell i Il965r) . we have derived two coupled differential equations 
governing the evolution of magnetic field perturbations in the presence of either ambipolar diffusion or the Hall effect. Solutions are found 
by WKB methods and by direct numerical integration. We find that instability arises from the combination of shear and non-ideal MHD 
processes, and is a result of the ability of these processes to open new pathways for the fluid to feed off the free energy of shear. They 
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(a) 



(b) 



Figure 8. Evolution of the initial state vector |ro) (solid arrow) under the transformation SR. (a) A counter-clockwise rotation by 6 is applied to |ro}. taking it 
into R|ro) (grey arrow), (b) This vector is then sheared along the j/-axis into \ri) (open arrow). The process represents the Hall-shear instabiliy (see text for 
details). Figure is drawn to scale for comparison with Fig.|2] 



tuin what would be simple linear-in-time growth due to current and vortex stretching from shear into exponential instabilities. We have also 
constructed a simple toy model based on transformation operators that not only captures all the qualitative results of this paper, but also 
matches the exact quantitative solution in some specific instances. 

In the case of ambipolar diffusion, anisotropic damping leads to the generation of magnetic field perturbations perpendicular to the 
background magnetic field. What ensues is a competition between ambipolar diffusion and shear, which is trying desperately to stretch the 
magnetic field along the stream-wise direction. In the end they both win, and the perturbations grow in a direction somewhere between the 
two that depends upon the geometry of the magnetic field and the ratio of the neutral-ion collision time-scale to the shearing time-scale. The 
resulting growth rates are on the order of 0.1 \2A\ (in the case of a time-independent background). It is notable that, in the general case of a 
time-dependent background, the growth rates increase in time without bound since highly-trailing (or -leading) shearing waves are trivially 
unstable. If the exponential growth phase lasts sufficiently long, growth rates may become comparable to or even exceed the shearing rate. 

In the case of the Hall effect, it is not current damping but rather current generation that gives rise to a magnetic field component 
perpendicular to the shear. Instability arises from the influence of shear on the propagation of circularly-polarized electromagnetic (whistler) 
waves, and is present so long as the shearing frequency \2A\ is larger than the ion cycl otron frequency (ti mes a factor proportional to the 
degree of ionization). Growth rates are on the order of 0.5 \2A\. In contrast to the work of I Wardld d 1 9991) and lSalbus & Tergueml ( 1200 II) . we 
find that instability depends not on ■ B or (fe • B) {k ■ CI), respectively, but rather on (fe • B) (fe • F), where F = V X v is the vorticity of 
the background. This should be negative for destabilization. Provided a given wavevector begins its evolution unstable (stable), it will remain 
unstable (stable) until non-linear processes intercede. The physical reason that typical growth rates for the Hall-shear instability are so much 
greater than the maximum growth rate for the ambipolar-diffusion-shear instability is that the Hall-shear instability employs a conservative 
process (cyclotron gyrations) rather than a dissipative process (ambipolar diffusion). The difference in these rates is related to the rate at 
which ion-neutral friction heats the gas. 

In both cases, unstable wavenumbers can be found for any sign of the velocity gradient. This e xplains why these processes were found to 
destabilize both inwardly- and outwardly-decreasing angular velocity gradients in accretion disks jBalbus & TerauenjboOlLlKunz & BalbusI 



20040 . While the analysis is complicated somewhat by the time-dependence of the shearing background, we find that the evolution unfolds 



as a series of time-independent problems, similar to studies of nonaxisymmetric instabilities in discs (e.g., Balbus & Hawley 1992b, Balbus 
& Terquem 2001). In the limit of small horizontal (y) wavenumber, ky/kz <C 1, the instability criterion may be written 

{k-VAf + k'^V-'^v <0, (58) 

with the resistivity tensor rj being given by either Equation l llOt (for ambipolar diffusion) or Equation l l37t (for the Hall effect). The final 
term here is the key to instability; in formal Cartesian index notation (i, j, k), it is k^rjijdvj /dxi. The maximum growth rate for a given rj is 
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, ^' ^ — . (59) 

For both non-ideal MHD effects, this is independent of the degree of ionization. Off-diagonal elements in the resistivity tensor are essential. 
No shear instabilities are present for isotropic damping processes, such as Ohmic dissipation. 

The impact of these results on astrophysical systems is unclear at this point. Our neglect of rotation precludes a straightforward appli- 
cation to accretion discs. There is a fundamental difference between planar shear flows and disc systems: in a planar shear flow there is only 
one characteristic gradient, dv jdx, whereas in a disk system there are two, one for the angular velocity, and one for the angular momentum. 
In astrophysical discs, the omitted Coriolis force generally dominates the shear dynamics. There is no asymptotic domain for either linear 
or non-linear perturbations in which the governing dynamical equations behave locally like Cartesian shear It seems that the best we can do 
here is to claim linear instability in constant-specific-angular-momentum discs. In addition, the well-known fact that inviscid planar shear 
layers are hosts to a plethora of non-linear instabilities suggests that the non-ideal MHD effects investigated in this paper may play only 
secondary roles. 

This may be undue pessimism, however. The instabilities investigated in this paper do have rotational counterparts, even in Keplerian 
systems. The real utility of this calculation, therefore, is that by removing rotation from the problem we obtain a clearer physical picture of 
what is going on in actual discs. One important simplification is that the MRI is absent. This is expected, of course, since the MRI is not just a 
shear instability. Rather, the MRI plays the crucial role of redistributing angular momentum and thereby opening paths to lower energy states 
in differentially-rotating discs. On the other hand, the ambipolar-diffusion- and Hall-shear instabilities have no preference for the source of 
the shear, whether it be in linear or angular velocity. The extent to which they lead to enhanced angular momentum transport is not governed 
by magnetic stresses, as in the case of the MRI, but rather whether or not the bulk fluid is responsive enough to utilize an increasingly radial 
magnetic field. 
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